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BEHAVIOR  OF  IONIZED  PLASMA  IN  THE 
HIGH  LATITUDE  TOPSIDE  IONOSPHERE 


I .  INTRODUCTION 

The  behavior  of  plasma  on  the  auroral  magnetic  field  lines  has  been 
the  subject  of  a  number  of  studies  in  recent  years.  Theoretical  models  of 
the  polar  wind  have  been  developed  by  Banks  and  Holzer  [1968,  1969],  Holzer 
et  al.,  [1971],  and  Lemaire  and  Scherer  [1973].  Using  the  13-moment  system 
of  transport  equations  of  Schunk  [1977],  Schunk  and  Watkins  [1981,  1982] 
have  studied  the  steady  state  flow  of  a  fully  ionized  H'*’  -  0*  -  electron 
plasma  along  geomagnetic  field  lines.  In  their  first  study  of  the  polar 
wind  [1981]  they  used  the  13-moment  equations  for  the  electron  gas  and  a 
simplified  set  of  transport  equations  for  the  hydrogen  and  oxygen  ions. 
Their  results  show  that  electron  temperature  anisotropy  develops  at 
altitudes  above  2500  km.  Below  2500  km  the  electron  gas  is  essentially 
collision  dominated.  In  their  next  study  of  the  polar  wind  [1982]  they 
have  used  the  13*'nioment  system  of  equations  for  the  hydrogen  ions  also. 
The  results  show  that  for  supersonic  flow  at  high  altitudes  the  hydrogen 


temperature  perpendicular  to  the  field  line  (T  ).  The  reverse  is  true  for 

Pi 

subsonic  flow. 

Using  the  same  13-moment  system  of  equations,  Mitchell  &  Palmadesso 

[1983]  developed  a  dynamic  numerical  model  of  the  plasma  along  an  auroral 

field  line.  The  plasma  consists  of  the  electrons,  hydrogen  and  oxygen 

ions.  The  electrons  and  the  hydrogen  ions  are  the  dynamic  species  in  the 

model.  They  have  performed  simulations  for  the  case  of  a  current-free 
MBnnaeript  apprtmd  May  23, 1986. 


polar  wind  and  the  case  in  which  an  upward  field  -  aligned  current  was 
applied  along  the  field  line.  The  results  of  the  polar  wind  simulations 
have  been  compared  to  those  of  Schunk  and  Watkins  [1981.  1982].  It  is  seen 
that  for  low  electron  temperature  and  supersonic  hydrogen  ion  outflow,  both 
exhibit  an  anisotropic  hydrogen  ion  cooling  and  both  have  a  region  in  which 
the  parallel  temperature  increases  with  altitude  before  adiabatic  cooling 
dominates.  The  electron  temperature  profiles,  however,  cannot  be  directly 
compared.  The  lower  boundary  of  Mitchell  and  Palmadesso  [1983]  was  held  at 
a  constant  temperature  and  the  simulation  was  run  until  a  steady-estate  was 
achieved.  The  outflow  condition  treats  the  upper  boundary  as  an  electron 
heat  sink  and  therefore  the  electron  temperature  decreases  with  altitude. 
Schunk  and  Watkins  [1981]  specified  the  electron  temperature  and  a  positive 
electron  temperature  gradient  at  the  lower  end,  as  the  initial  condition  to 
study  the  electron  temperature  anisotropy.  The  upper  boundary  need  not  be 
specified  in  this  case. 

We  have  developed  a  numerical  model  to  study  the  steady  state  behavior 
of  the  plasma  encompassing  geomagnetic  field  lines.  The  theoretical 
formulation  used  is  based  on  the  l6~moment  system  of  transport  equations  of 
Schunk  et  al.  [1982].  The  13*^oment  system  of  transport  equations  allow 
for  different  species  temperature  parallel  and  perpendicular  to  the  field 
line,  but  allow  only  a  single  heat  flow  per  species.  The  16-moment 
equations  allow  transverse  and  parallel  thermal  energy  to  be  transported 
separately,  which  simulates  the  behavior  of  a  large  temperature  anisotropy 
for  a  collisionless  plasma  better  than  the  13‘*’moment  equations.  This  is 
expected  to  be  an  Important  advantage  in  studies  of  auroral  field  aligned 
transport  which  include  the  effects  of  return  currents  and  anomalous 
transport  processes,  such  as  anomalous  resistivity  and  associated 


anisotropic  heating.  Such  studies  are  now  in  progress,  and  will  be  the 
subject  of  a  future  report.  We  have  used  this  model  to  perform  polar  wind 
simulations.  This  to  our  knowledge  is  the  first  successful  steady  state 
solution  to  the  16  moment  set  of  transport  equations  for  the  polar  wind 
problem . 

II.  THE  MODEL 

A  numerical  model  has  been  developed  to  simulate  the  steady  state 
behavior  of  a  fully  ionized  plasma  (H"*",  O’*”  and  the  electrons)  along  the 
geomagnetic  field  lines  in  the  high  latitude  topside  ionosphere.  The 
electrons  and  the  hydrogen  ions  are  the  dynamic  species  in  the  model  and 
the  oxygen  ions  form  a  static  background  population  at  a  constant 
temperature.  The  H*  and  electron  temperature  anisotropies  and  heat  flows 
are  not  sensitive  to  the  oxygen  temperature  and  heat  flows,  as  noted  by 
Schunk  and  Watkins  [1981,  1982].  Thus  maintaining  the  0*^  ions  at  a 

constant  temperature  did  not  Introduce  appreciable  errors  as  far  as  the 
ions  and  the  electrons  are  conerned  [Schunk  1982]. 

The  distribution  function  is  considered  to  be  gyrotropic  about  the 
field  line  direction,  which  reduces  the  1 6-moments  approximation  to  six 
moments:  number  density,  velocity  of  the  species  parallel  to  the  field 
line,  temperature  pa.  allel  and  perpendicular  to  the  field  line,  and 
parallel  and  perpendicular  heat  flow  along  the  field  line.  The  resulting 
transport  equations  used  for  hydrogen  and  electrons  are  given  as  follows 
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For  a  given  moment  F  of  the  distribution  function,  6F/6t  represents 
the  change  in  F  due  to  the  effects  of  collisions  and  may  also  include 
anomalous  transport  effects  associated  with  plasma  turbulence.  The 
collision  terms  used  in  this  model  are  those  of  Mitchell  and  Palmadesso 
[1983] •  They  have  used  Burger's  [1979]  collision  terms  for  the  case  of 
Coulomb  collisions  with  corrections  for  finite  species  velocity 
differences.  The  collision  terms  are  given  in  Appendix  A. 

We  have  assumed  that  the  total  flux  tube  current  I  remains  constant 
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which  implies 
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In  the  following,  the  current  I  Is  assumed  to  be  zero 


Using  equations  (2),  (7)  and  (8)  the  electric  field  E  parallel  to  the 
field  line  is  calculated  [Mitchell  and  Palmadesso,  1983], 
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We  will  consider  the  steady  state  solutions  of  the  equations  (1)" 
(6).  This  set  of  ordinary  differential  equations  are  solved  by  first 
eliminating  all  derivatives  In  a  given  equation  except  one  and  then 
Integrating  the  resulting  quatltles  numerically  along  the  field  line.  The 
new  set  of  equations  are: 
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The  simulation  Is  carried  cut  on  an  unequally  spaced  grid,  with 
smaller  cell  sizes  at  the  lower  end  of  the  flux  tube.  This  enables  us  to 
determine  the  plasma  dynamics  properly  in  the  presence  of  large  density 
gradients  due  to  the  small  scale  height  of  oxygen. 


The  equations  we  have  to  solve  represents  a  system  of  stiff 
differential  equations.  A  typical  solution  of  a  stiff  system  has  a  short 
initial  interval  in  vrtilch  it  changes  rapidly  (called  an  initial  transient) 
after  which  it  settles  down  to  a  comparatively  slowly  varying  state. 
Methods  of  solving  stiff  ordinary  differential  equations  has  been  discussed 
by  Gear  et  al.,  [1979].  We  will  use  the  general  procedure  used  by  Watkins 
[1981]. 


III.  SIMULATIONS  AND  RESULTS 


Schunk  and  Watkins  [1982]  have  performed  polar  wind  simulations  for 
both  supersonic  and  subsonic  hydrogen  ion  outflow.  Until  recently  it  was 
not  clear  if  the  polar  wind  la  supersonic  or  subsonic.  Nagai  et  al.  [198*1] 
first  reported  the  supersonic  nature  of  the  polar  wind  along  polar  cap 
field  lines.  The  data  for  this  study  were  obtained  from  DE  1  satellite. 


This  pap'>r  clearly  shows  that  the  hydrogen  ion  flow  in  the  polar  ionosphere 
is  supersonic  in  nature.  Subsonic  polar  wind  has  not  yet  been  observed. 
We  shall  thus  consider  hydrogen  ion  outflow  as  supersonic 


As  noted  above,  solutions  have  been  obtained  for  the  current  free 
case.  The  lower  boundary  is  fixed  at  1500  km.  At  this  boundary  we  set  H'*’ 
ion  velocity  at  16  km  3“^,  H'*’  ion  density  at  80  cm'^,  and  oxygen  ion 
density  at  5000  cm”^.  The  hydrogen  ion  temperature  at  the  lower  boundary 
is  3500®K.  The  oxygen  ion  temperature  was  kept  constant  at  1200®K,  The 
results  of  this  current  free  polar  wind  simulation  are  shown  in  figures  1a- 
1g.  Figure  la  shows  oxygen  ion  density,  1b  shows  H"*"  ion  density,  Ic  shows 
electron  density.  Id  shows  hydrogen  ion  velocity.  Oxygen  ions  are  the 
dominant  species  only  up  to  an  altitude  of  3500  km.  Supersonic  H'*'  ions 
flow  in  divergent  magnetic  field  lines.  The  hydrogen  momentum  balance 
equation  shows  that  the  amblpolar  electric  field  value  produces  the  sharp 
increase  in  ion  velocity  at  the  lower  end  of  the  tube.  Electron 


velocity  also  Increases  as  is  evident  from  equation  (9)  and  is  shown  in 
Figure  1e.  The  total  flux  along  the  flux  tube  is  conserved,  i.e.,  nVA  - 
constant.  Therefore,  as  the  velocity  and  the  area  Increases  the  density 
should  decrease.  The  hydrogen  ion  density  decreases  as  it  flows  through 
the  diverging  flux  tube.  The  electron  density  also  decreases,  as  is 
evident  from  the  charge  neutrality  equation  (7). 

A  study  of  equation  (11)  reveals  that  at  the  upper  end  of  the  flux 
tube  the  velocity  variation  is  driven  mainly  by  the  terms  containing  the 


hydrogen  ion  temperature  difference  (Tp^  -  electron  temperature 

difference  (T  -  T  ),  hydrogen  ion  temperature  (T^  )  and  electron 
®l  ®i  P| 

temperature  (T  )  parallel  to  the  field  line,  i.e.,  mirror  forces  and 

®  I 

parallel  pressure  gradients. 


Figure  If  shows  the  hydrogen  ion  temperature  parallel  and 
perpendicular  to  the  geomagnetic  field.  The  hydrogen  Ion  temperature 
exhibits  two  basic  characteristics: 

(1)  Adiabatic  cooling  -  Supersonic  Ion  gas  cools  down  as  It  expands  In  a 
diverging  magnetic  field 

(2)  Temperature  anisotropy 

These  effects  were  also  displayed  by  the  previous  studies  of  the  polar 
wind  using  the  13-moment  system  of  equations  by  Schunk  and  Watkins  [1982] 
and  Mitchell  and  Palmadesso  [1982].  Like  Schunk  and  Watkins  [1982],  we 
obtained  solutions  for  only  positive  heat  flow  at  the  lower  boundary. 
This  Implies  that  an  upward  flow  of  heat  from  the  lower  Ionosphere  Is 
required  for  a  supersonic  hydrogen  Ion  outflow  [Schunk,  1982].  This  upward 
flow  of  heat  Is  associated  with  a  negative  temperature  gradient. 

At  the  lower  end  of  the  flux  tube,  due  to  rapid  expansion  of  the 
supersonic  Ion  gas,  both  T^  and  T^  decrease.  This  was  also  noted  by  Schunk 
and  Watkins  [1982].  The  hydrogen  Ion  temperature  perpendicular  to  the 
field  line  decreases  continuously.  Appreciable  temperature  anisotropy 
develops  around  2500  km  .  The  temperature  anisotropy  Is  more  prominent  In 
the  Ions  than  In  the  electrons.  This  Is  because  of  the  fact  that  the 
hydrogen  Ions  have  supersonic  velocities  and  In  the  colllslonless  region 
heat  transport  Is  not  efficient  for  Ions,  compared  to  that  for 
electrons.  Also,  as  the  particles  move  up  In  a  decreasing  magnetic 


field,  V  decreases  In  order  to  keep  adiabatic  Invariant  p  constant.  As 

1  2 

the  total  energy  of  the  particle  Is  conserved,  ^  mv  decreases  with  the 
1  2 

Increase  of  2  energy  gets  converted  to 

parallel  energy.  A  study  of  equation  (111)  and  (Iv)  show  that  the 


temperature  variation  Is  produced  mainly  as  a  balance  between  the 


convection  process  and  the  mirror  force  effects.  The  collisions  play  an 
Important  role  below  2500  km.  In  equations  (v)  and  (vl)  the  temperature 
gradient  terms  are  very  important. 

We  compared  our  solution  of  the  equations  based  on  the  1 6-moment 
system  with  that  of  Schunk  and  Watkins  [1982]  solution  of  the  13~nioment 
equations.  The  general  behavior  of  the  ion  temperature  parallel  and 
perpendicular  to  the  field  line  remains  the  same.  Ion  temperature 

perpendicular  to  the  field  line  decreases  continuously  with  altitude. 
Parallel  temperature  first  decreases  and  then  Increases  and  finally  tends 
to  remain  constant  with  altitude.  We  observe  an  anisotropy 

with  T^  >  T|  below  2500  km.  Schunk  and  Watkins  [1982]  observed  similar 
anisotropy  for  the  low  electron  temperature  case  when  they  lowered  the  H* 
Ion  temperature  and  heat  flow  and  when  the  oxygen  Ion  density  was  increased 
at  the  lower  end.  They  have  pointed  out  that  the  results  are  sensitive  to 
the  boundary  conditions  used,  in  the  sense  that  the  qualitative  picture 
remains  the  same  only  the  quantitative  picture  differs. 

We  have  also  compared  our  results  with  those  of  Mitchell  and 
Palmadesso  [1983]  and  Holzer  et  al.  [1971].  Mitchell  and  Palmadesso  [1983] 
started  their  simulation  at  800  km  at  the  lower  boundary.  At  this  level 
the  H'*’  ion  temperature  Increases  with  altitude  due  to  Joule  heating  caused 
by  ion  collisions  with  oxygen  ions.  This  effect  rapidly  decreases  with 
altitude.  The  general  qualitative  picture  here  is  in  good  agreement  with 
our  results. 

Holzer  et  al.  [19tl]  have  compared  the  solutions  of  hydrodynamic  and 
kinetic  solutions  for  supersonic  polar  wind  outflow  in  a  collisionless 
regime.  Both  hydrodynamic  and  kinetic  solutions  show  that  T^  decreases 
continuously.  T^  for  both  cases  decreases  rapidly  with  altitude  and  then 
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tends  to  become  constant.  For  the  hydrodynamic  case  at  all  altitudes  the 
hydrogen  Ion  temperature  parallel  to  the  field  line  was  greater  than  that 
perpendicular  to  the  field  line.  The  kinetic  solutions  showed  a  region  at 


the  lower  end  where  was  greater  than  .  The  cross-over  region  was 
noticed  at  the  lower  end.  In  these  cases  also  we  notice  that  our  solutions 
are  In  good  agreement  with  Holzer  et  al,  [1971]. 

Finally,  we  compared  our  results  with  Nagal  et  al.  [1984].  Nagal  et 
al.  [1981]  have  used  the  data  from  Dynamics  Explorer  (DE)  to  demonstrate 
the  supersonic  nature  of  the  polar  wind.  The  observations  were  obtained 
from  65**  to  8l**  Invariant  latltute  and  an  altitude  near  2R£.  The  results 
show  that  for  an  estimated  range  of  spacecraft  potential  of  3  to  -^5V,  the 
temperature  range  of  0.1  to  0.2  eV  was  obtained  corresponding  to  flow 
velocities  of  25  to  16  km  s'^  Calculated  Mach  numbers  ranged  from  5.1  to 
2.6.  The  Mach  number  calculated  from  our  results  corresponding  to  a  flow 
velocity  of  16  km  per  second  Is  2.2. 

We  have  varied  slightly  the  ion  velocity,  density  and  temperature 
at  the  lower  end  to  study  the  sensitivity  of  the  solutions,  to  these 
parameters.  At  the  lower  end  the  hydrogen  Ion  density  Is  Increased  to  100 
cm~^,  velocity  Is  decreased  to  10  km  per  second  and  the  temperature  Is 
raised  to  3800**K.  The  results  of  these  changes  are  presented  in  figures 
2a-2f .  The  qualitative  picture  and  the  relative  Importance  of  the 
contributing  terms  remain  the  same  as  before.  The  quantitative  change  Is 
obvious  as  the  Initial  conditions  are  changed.  The  mach  number  calculated 
in  this  case  of  a  flow  velocity  of  16  km  s"^  is  2.68.  This  is  In  good 
agreement  with  the  value  2.6  calculated  by  Nagal  et  al.  [1984]. 

Electron  temperature  parallel  and  perpendicular  to  the  geomagnetic 
field  lines  are  shown  In  Figure  1g.  The  most  Important  effect  exhibited  by 


electron  temperature  la  : 

( 1 }  temperature  anisotropy 

(2)  Below  2500  km  the  electron  gas  is  collision  dominated,  as  noted  by 
Schunk  et  al.  [1981] 

"Collision  dominated"  Implies  that  the  collision  relaxation  rate  Is  greater 
than  the  rate  associated  with  the  driving  term  that  produces  the 
anisotropy.  Appreciable  temperature  anisotropy  develops  above  2500  km. 
Both  temperature  parallel  and  perpendicular  to  the  field  line  Increases. 
Heat  flow  Is  considered  downward  here  (that  Is,  heat  flows  from  the 
magnetosphere  to  the  Ionosphere)  as  considered  by  Schunk  et  al.  [1981]. 
Also  as  noticed  by  Schunk  et  al.  [1981],  the  electron  temperature  ratio  at 

the  upper  end  la  T  /T  -  2  .  The  contributions  from  the  collision  term 

i  I 

are  very  Important  below  2500  km  In  equations  (111)  and  (Iv).  -Since  the 

collision  terms  are  Important  In  this  region,  the  electron  temperature 

anisotropy  observed  In  this  region  Is  very  small.  The  anisotropy  Increases 

with  altitude.  Electron  temperature  variation  Is  produced  mainly  as  a 

balance  between  the  convection  process  and  the  mirror  force  process.  The 

heat  flow  Is  downward  for  the  electrons  and  the  perpendicular  energy 

Increases  as  the  parallel  energy  decreases.  The  temperature  gradient  terms 

are  Important  at  all  altitudes  In  equations  (v)  and  (vl)  the  electron 

species  velocity  Is  less  than  the  electron  thermal  velocity. 

Our  results  here  are  In  good  agreement  with  Schunk  and  Watkins  [1981] 

low  electron  temperature  case.  Schunk  and  Watkins  [1981]  varied  the 

« 

boundary  electron  temperature  gradients  to  see  the  extent  to  which  their 
solutions  are  valid.  Their  results  show  that  with  the  Increase  In  boundary 
electron  temperature  gradient,  the  electron  temperature  Increased  but  the 


direction  of  anisotropy  remained  the  same.  It  Is  also  clearly  explained  by 
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Schunk  and  Watkins  [1981]  that  this  observed  anisotropy  is  opposite  to  that 
predicted  by  the  Chew-Goldgerger-Low  (  double  adiabatic  )  energy  equations. 
The  tendency  of  the  adiabatic  terms  to  produce  an  anisotropy  with 
is  also  present  in  our  equations,  but  the  heat  flow  effects  dominate  and 
act  to  produce  the  reverse  anisotropy. 

The  electron  temperature  profiles  cannot  be  directly  compared  to 
Mitchell  and  Palmadesso  [1983].  They  treated  the  upper  boundary  as  an 
electron  heat  sink,  which  decreases  the  electron  temperature  with  altitude. 


IV.  DISCUSSIONS 


We  have  studied  the  steady  state  behavior  of  the  plasma  encompassing 
the  geomagnetic  field  lines  using  equations  based  on  the  1 6-moment  system 
of  transport  equations.  The  asymptotic  behavior  of  the  solutions  of  the 
transport  equations  at  high  altitudes  is  quite  sensitive  to  the  boundary 
conditions  applied  at  the  lower  end  of  the  flux  tube.  This,  coupled  with 
with  absence  of  any  standard  set  of  assumptions  for  specifying  those 
dynamic  variables  for-  which  upper  boundary  conditions  must  be  supplied, 
leads  to  considerable  variation  of  the  results  at  high  altitudes  in  the 
studies  referenced  above.  Effectively,  one  finds  similar  ionospheres  and 
widely  different  magnetospheres  in  this  group  of  calculations.  While  this 
makes  quantitative  comparison  of  results  somewhat  difficult,  it  is  clear 
that  qualitative  agreement  among  the  various  calculations  is  good. 

Studies  of  auroral  field  line  equilibria  including  the  effects  of 
return  currents  and  anomalous  transport  processes,  and  using  the  16  moment 
approach  described  here,  are  in  progress. 
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Appendix  1 
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